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ABSTEACT 



The object of this report is to investigate the response of a shaft 
or beam to shock applied transversely at the end supports. As a secondary 
objective, a simplified method of designing a shaft for shock resistance 
is suggested. In particular, the results of this work demonstrate the 
prediction of the transient history, with respect to time and position, 
of the beam deflection, shear force, and bending moment and show a calciu- 
lation of the times, positions, and maxirnum amplitudes of these quantities. 

To accomplish the objective, velocity shock was assumed and the 
differential equation for flexural vibration of Timoshenko was solved. 

The solution was then placed in a form such that values depend only on 

the ratio of the beam length to radius of gyration of a cross section and 

the ratio of shear modulus, modified by a factor accounting for section 
shape, to the modulus of elasticity. The maxima for deflection, shear and 
bending moment were then tabulated for various length to radius of gyra- 
tion ratios. 

As a model, only a solid, circular, steel shaft of six inch radius 
was used, but the method of solution could be extended to other config- 
urations . 

As a result of the investigation, it can be shown that, comparing all 
simply supported uniform beams of constant ratios of length to radius of 
gyration of a cross section and of constant ratio of shear modulus to 
Young’s modulus, modified by section shape, the time histories of deflec- 
tion, shear, and bending moment will differ only in magnitudes, i.e., the 

histories of any such beam can be scaled from that of another. 

It was found that there is a range of length to radius of gyration 
ratios where rotary inertia of elemental cross sections and shear effects 
become significant. Further, there are a two fold infinity of natural 
frequencies, each of which corresponds to an odd numbered deflection mode. 
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The higher series of frequencies is associated with axial compression 
and tension due to bending^ and the lower series is associated with the 
shearing effect* It was found that the higher frequency series^ for all 
modes, goes to infinity as the length to radius of gyration ratio goes 
to infinity and that, as a resiolt, the equations derived reduce to the 
classical flexural vibration equations. 

It is believed that the results obtained in this work coiiLd provide 
a designer concerned about the effects of shock with a useful tool in 
anticipating design performance due to shock. 

There was an indication that the energy contained in a beam subjected 
to shock cannot be adequately described with considerations of overall 
deflection. However, this point should be investigated further. 

It is recommended that the equations and curves obtained be checked 
by experiment. It was noted that the equation for shear converged rather 
slowly; so, even though the curves obtained from coir^uter calculations 
included terms through the forty-ninth mode, this may not provide enough 
accuracy. It is further recommended that other types of supports, 
especially fixed ends and cantilever, be investigated. 
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A = cross sectional area of beam. (in. ) 

E = modulus of elasticity, (psi) 

G' = Shear modiolus, (psi) 

G = k'G' (psi) 

2 

g = acceleration due to gravity (in. /sec. ) 

I - moment of inertia of a cross section about the neutral axis, (in.^) 

J = mass moment of inertia of a unit length of shaft about the neutral 
axis. (lb. -in. -sec. 2/in. ) 

k’ = factor relating average shear stress to maximum shear stress. 



, rnr /. -1\ 

k = = — . (in. ) 

L = length of beam. (in. ) 

M = bending moment, (in. -lbs.) 

2 2 

m = mass per unit length (lb. -sec. /in. ) 
n = mode number, an integer (n = 1,3^5>****) 

Q = shear force, (lbs.) 
r = radius, (in. ) 
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= shear stress, (psi) 

Tin^ T^n = time period of the nth mode of the lower and the higher series 
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V = initial velocity, (in. /sec.) 

tUin^ UJ^n ~ frequencies of the nth mode corresponding to the lower and the 
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X = position along shaft, measured from one end. (in. ) 



V 



y = transverse deflection. 

Q 

z = dummy variable = 1 + — 

JIj 

2 

>, j /•lb. -sec. / . 3\ 

p = density ( / in. ) . 




Subscripts : 

n = mode number = 1,3;5;*«* 
b = refers to bending only, 
s = refers to shear only. 

Primes, except in k’ and r', indicate dimensionless values. 
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I. ITJTRODUCTION 



Today ^ through use of such devices as armor plating^ subdivision, 
and torpedo bulkhead systems, naval ships are well able to remain afloat 
after being struck by sizeable amounts of explosives. However, a ship is 
useless if the shock associated with an explosion causes vital machinery 
to become inoperable. Furthermore, with modem nuclear weapons, depth 
charges, and mines available, physical contact with the ship is not 
always possible; and shock effects are then the primary means of incapa- 
citating ships. This is particilLarly true for submarines. Indeed in any 
industrial plant or ship where shock or jolting is anticipated, a design 
must be able to tolerate these effects in order to perform its mission. 

Of particular susceptability to the effects of shock are shafts and 
beams when their bearings or supports accelerate perpendicular to the 
shaft axis. In actual situations, such items as piping, condense»r shells, 
rotating shafts, and in some cases plates could fall into the general 
category of beams. Of interest to the engineer are such problems as the 
amount of deflection that must be anticipated in order to have sufficient 
clearance around the shaft, and the magnitude of the stresses involved 
in order to assure that no plastic deformation occurs. Thus, this inves- 
tigation attempts to predict these values. 

There have been many works on the dynamic flexural behavior of 
beams. Among the primary of these is that of Timoshenko [l]. However, 
these works neglect the effects of shear and rotary inertia which becomes 
significant when the length to radius of gyration ratio becomes small. 
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To the author knowledge, there is no simple design method established 
to account for shock effects. 

The investigation was limited to a solid, circular beam, simply supported. 
The effects of spinning, such as in rotating machinery, on shaft response was 
considered negligible; and no calculations account for this. 

The excitation was considered to be a "velocity shock", i.e. the beam 
supports experience an instantaneous constant velocity at the incidence of 
the shock. This is a commonly used way of specifying shock excitation. Use 
of velocity shock allows calculations to be simplified; for when a body, 
initially at rest, is excited by a velocity shock, its response relative to 
the point of excitation is precisely the same as if the body were initially 
in motion with a velocity equal to the excitation velocity, and suddenly 
stopped m . Thus calculations can be made as though the beam supports are 
stationary and an initial constant velocity exists along the length of the 
beam. 

If one looks at a beam subjected to shock in actual practice, the response 

woiild depend upon the direction from which the excitation eminated. Thus, if 

one support were excited before the other, the response would be different 

than if both ends were excited at the same time. In an atteri5)t to simplify 

difficult calcinations, this report is limited to the case of symmetrical 

*• 

excitation. 

Linearity has been assumed in utilizing the equations involved in this 
report. Thus, the results are limited to beams whose stresses always remain 
in the elastic region. It was felt that this would not place a serious 
handicap on a designer if it is assumed that a beam has failed when the stresses 
resulting from calculations demonstrated herein exceed the yield limit of the 
material involved. 
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To Siam up, then, the model used for the investigation is a solid circular 
beam, simply supported, with ends fixed; it is subjected to a velocity shock 
along its entire length. The response of the model is limited to the elastic 
limits of the beam material, which was assumed to be steel. The object of 
the work is to describe the response of the beam model with respect to deflec- 
tion, shear, and bending moment; then to show to what limits the results apply 
to other beams than the model chosen; and lastly to propose a simplified method 
of calculating the deflection and stresses for other beams. 
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II. EROCEDURE 



Definition of Problem . Whereas the physical situation specifies that, 
at time, t = o, the ends attain an instantaneous velocity, v, the problem 
shall be considered as one where, at t = o, the entire length of the beam 
attains a velocity, -v. The final equation for deflection, then, will 
be the sum of the values derived eind vt. 

Figure I indicates the model investigated and the sign conventions. 
Motion is restricted to the x,y plane. The problem is to find the motion 
of the beam Indicated in Figure I. 



Derivation of the Differential Equation of Motion^ Assume that the 
deflection is due to two separate and separable effects, shear and bending. 
Then total deflection is 



y = Yg + , (1) 

where 

y = deflection due to shear only 
s 

y^ = deflection due to bending only. 



Look now at Figure I. The elemental mass is considered rigid. 
Consider that, for all shear deflections, there is no rotation of the 
element. Rotation is associated with bending deflection only. Thus 






V oM •^b 

/. Moments = Qdx - dx + J — « — dx = 0 , 



where the last term represents rotary inertia of the elemental mass. Thus 
it can be seen that, for bending moments, the element always remains 



Figure I. 

Beam Model and Sign Convention 
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M + dx 






'f 

mgdx 
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perpendicular to the slope of the deflection due to bending, ; 

ox 

— ^ simply represents the angular acceleration of the element. 
3t Sx 



From this it is clear that 






9 t 9x 



Equating forces in the y direction, 



I 



90 

forces = - - 5 ^ dx 

9x 



mdx ^ = 0 , 

9t 



or 



^ - m ^ 

9t" 



Recalling that the maximum shear stress. 



^s(max) 



= G’i , 



where \|i represents the slope of deflection due to shear, then 



. ^^s 

S , ,. = -G — 

s(inax) ox 



Although shear stress is not constant over an entire cross section. 



be shown that the average shear stress, S , is 

s 

S = - k G ^ = - G ^ , 
s ox 9x 



and 



( 2 ) 



( 3 ) 



it can 
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where k is a factor depending upon section shape H • 



Thus 



Q = S^ A 



= - GA 



5x 



Continuing, it is well known that 



M = + El 



9x^ 



(^) 



(5) 



Combining equations (l), (2), (3), (^), and (5); one obtains 



El 






, 1 * 

Elm , ° ^ 



+ m 






Jm 

GA 



aV 



= 0 



( 6 ) 



Equation (6) is the equation to be solved. Its derivation is 
presented in reference H , but has been given here in slightly different 
form because use will be made of the equations Involved in the derivation. 



Solution of Equation of Motion . Assume that the variables are separable. 
Then 

y (x,t) = T (t) X (x) . 

The equation in x is then of the form 

CX. 0, 

dx dx ^ 

where A, B, and C are functions of time. The solution of equations of this 
type is of the form 

X (x) = C^ sin k X + Cg cos k x + C^ sin h k x + Cij^ cos h k x . (7) 
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Since the ends are pinned, deflection and bending moments are always 
zero at the ends . Thus the boundary conditions are 



X(0) = Xg(0) + X^(0) = 0 



X(L) = Xg(L) + X^(L) = 



(8) 






dx 



x:=o 






dx 



= 0 



x;=L 



Putting equations (8) into (7) it can be shown that X (x) is of the 
form, 



X (x) = 



GO 

z 



n=l,3,5... 



C Sin k X , 
n n 



(9) 



where k 

n 



ntT 

L 



n = 1,3,5, 



See Appendix A for proof of this. 

Assuming no damping, a solution for T (t) can be of the form, 

T (t) = C, cos to t + C„ Sin to t , 

n ^ In n 2n n 

or 

y(x,t) 

Putting (lO) into (6) it is easily seen that the frequency equation results, 
i.e. , 



00 

r [0 

n— 1 ,3,5,»** I — 



In 



cos (0 t + Sin to t sin k 



n 



2n 



n 



] 



n 



( 10 ) 




(11) 
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Equation (ll) shows that there are two frequencies associated with 
each which fact will be amplified in a later section. Thus, the 
assumed T(t) is incomplete, since it does not account for both frequencies. 
Allowing for both frequencies, let 



y(x,t) = Y1 

n=l^ 3^ 5 y ' 






- sin UL t + C_ cos ti)_ t + C- sin t 
In in 2n In 3n 2n 



^ “an* 



; sin k 

J. n 



X . 



( 12 ) 



From equation (l2), it is seen that four initial conditions must be 
utilized to evaluate the constants. Thus, initially, the deflection of 
the beam is zero. From the definition of the problem, the initial 
velocity along the entire beam is the velocity of shock, -v . Lastly, 
due to rotary inertia, the slope and angular velocity of the deflection 
caused by bending cannot instantaneously change, and are thus zero • 
Putting this in mathematical terms. 



y (x,0) = 0 



dt 



t=o 



“V 



It-O 



= 0 



dxdt 



= 0 



t=o 



(13) 

(14) 



(15) 



(16) 
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Now, put equation (12) in (3) and integrate. Then 



Q = 



dQ 



dx = -m 



iz 



dx = - 



Z 

n— 1 ^3^5^ ••• 



^ I C, sin W, t 

kn in In In 



2 2 2 ~1 
+ w. cos CO. t + co„ C_ sin ou t + uu Ci, cos u) ^ cos k x 

In 2n In 2n 3n ^n ^n 4n 2n I n 



+ Const. 



(17) 



But from Appendix A, it can be seen that, at x = ^, Q = 0, thus const. = 0. 
Apply equation (4) to (iZ), so that 



^s = 



r , 1 

= ■ GA 



Qix = 



i i'B 



n=l,3^... n 



‘^in ®2n ’-"in* “2n =3n “2n’^ * “2n* I V' 






(18) 



Again, boundary conditions preclude constants of integration. Thus5 from 
equations (l):, (l2), and (l8). 



^b = 



n— 3^ 5.? • • • 



irf^ 



1 - 



In 



GAk. 



n 



sin to; t + cos (0-, t 
In In 2n In 



2 




— . — 


\ 


, “^2n 








^ “ 2 

GAk 

- n_ 




°3n “an* '^l.n “2n‘ 











sin Ic^x. (19) 



Now the initial conditions can be applied. Applying equation (l3) and 
looking at the nth term. 
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(20) 



(x.O) = (Cgn + V " 



^2n = -^4n . 



Applying equation (l4) , 






n 



at 



t=o 



' (®ln “in * ^3n ®2n> >« - -V , 



or (C. u). + C_ tu„ ) 

' In In 3n 2n^^ 






(-v) sin k xdx = - — 
^ n rnr 



( 21 ) 



Continuing, from equation (15 ) ^ 



9x 



= k 



t=o 



n 



muo. 



^2n 



In 



GAk 



n 



ittu. 



" ^4n '"■ 



2n 



GAk 



n 



cos k x= 0. 
n 



But from equation (20), this means that, 



nW- rai)^ 

In 2n 



"2n 



= 0 



GAk GAk 

n n 



Since, as will be seen later in this work, cannot equal ^ 

contradiction results, and 



^2n = -^n = 



( 22 ) 



Lastly, apply equation (l6), so that. 



3^yi, 



axB| 



= k 



n 



t — o 



mP. 



\n^n 



In 






GAk 

n 



^ ^3n '«2n* ^ - 



2n 



GAk 



n 



cos k X 
n 



= 0 . 



(23) 
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By combining equations (21 ) and (23), 



C 



In 



W 

nn 



(GAk ^ - mu)_ 

' n 2n ' 

( 2 2v 

mco^ (co„ - J 
In ' 2n In'' 



C 



3n 



nTT 



(GAk 2 
n 




rm, 



2n 





(24) 



(25) 



Thus, by applying equations (22), (24) and (25) to equation (l2) and 
adding the quantity, vt, as indicated in Definition of Problem above, 
one finally obtains 



y(x,t) 





n— 1 , 3 ^ 5 • • • 



(GAk 2 muj, ) 

' n - m In ^ 

- 2 j 

(cu- - ) 

2n 2n In'^ 



sin UJ- t 
2n 



(GAk ^ - ma)_ ^) 

' n 2n ' 

• ^ 2 
moj^ (u)„ - ) 

In ' 2n In-^ 



sin 



V 



sin k X 
n 



(26) 



Solution of the Frequency Equation . To complete the solution of the 

differential equation of motion^ the tvo series of frequencies must be 

evaluated^ i^e. equation (ll) must be satisfied. 

2 

Thus solving equation (ll) for o)^ , one finds that 




Clearly, there are two frequencies associated with each k^ ; thus, there 
is a two-fold infinity of frequencies, each pair of which is associated 
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with an odd mombered mode of deflection. Intuitively^ since the deflection 
was broken into components from shear and from bending^ one would expect 
one of these frequency series to arise from the bending and the other 
from the shear. This point will be demonstrated more rigorously in a 
later section. 

It has been found more convenient to put equation (27) into the form 




(2Ta) 



or, letting 




2GAk ^ 



(2Tb) 



Details of this derivation are given in Appendix B. 



Calculation of Shear and Bending Moment . Based on the equations derived 
in the previous sections^ the values of shear and bending moment may now 
be calculated. 

Combining equations (22)^ (24), (25)^ and (l7) the shear force at 
any section, at any time, is 
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Q = 



y~ !n 

L — riTT 
n=l ^3^5* • • 



(GAk/ - 



k («)/ - ID, I 
n ^ 2n In ^ 



?s '"2n “2n’= 



2 2 
( GAk - imu„ 

n 2n . + 

7! 7 ? ou, sin (ju, t 

2 X In In 



k ((Ho - ) 

n ' 2n In ' 



cos k X 
n 



( 28 ) 



Utilizing equations (2), (19) and (28), it can be shown that, 



5hf _ V W 

Sx ~ /_ nrr 



3^ 5 • • ' I 



""2n „ ™2n 

* - GAk 



n 



n 



GAk - mou, 
n In 

r 2 2 

m (uJo - 00 , 

' 2n In 



sin (Ho t 
2n 



imH. 



In 



moo. 



+ J k. 00, 



In 



k I n In GAk 

n \ n 



GAk - moOo 
n ^ 

“ <"2^ - 



Sin cu, t > cos k x . 
In ( n 



Integrating this. 



M = 



W 



nrrk 



n-1,3,5- 



n 



moo, 



2n 



mcH, 



+ J I k 00^ 



2n 



k ' " \ " n 2n GAk 
n \ n 



2 2 
GAk - mo), 
n In 

2 57 

m (cOo - 00 , ) 
^ 2n In 



sin oOo t 
2n 



m£B. 



In 



m(U. 



+ J I k 00, 



In 



k \ n In GAk 

n \ n 






GAk - mcu, 
n In 

2 2 

m ((Ho - 00 , ) 
' 2n In 



sin oo^t V sin k^x , ( 29 ) 



where once again boixndry conditions preclude a constant of integration. 



Evaluation of Equations . An indication of the validity of any equation 
occurs when the equation reduces to a simpler, classically known form 
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based on restrictions of the parameters. Accordingly, look at the values 
of y, Q, and M when L/r' goes to infinity. 

Thus, from equation (SJa), 



lim 

L 



(U = 
n 



m 



2GAk,- 



G /L 

E („ > 

Vrrr n 



^ (rrr'n ) 



from which it can be seen that. 






2n 



-> “ 



(U. 



ln“ 



T.A, 2 2 2 .2. 

EAk n n- r 
n 



mL 



k** Si 

n m 



or 






■2n—> ’ “iS ? ■'n J— ■ 



(30 a,b) 



Notice that is the classical natural frequency for beams of long 
length to radius ratio. Putting equations (30 a,b) into (26) then 




But this is just the classical deflection equation for large ~i |]lj 

Now, look at the bending moment equation, equation (29). At first 
glance, it appears that the coefficient of the sin term blows up 
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as ^ — > ® . However, it can be shown that this term actually goes to 
zero. This is because the mass moment of inertia, J, times goes to 
zero. See Appendix C for proof of this convergence. 

Appendix C also indicates that 



lim 



M 



r~ kv “‘“in . . . V 

= > — — Sin t Sin k x , 

£ n 2 In n ' 



— j — > “ n=l , 3 > 5 • • • 

r 



n 



or, for large — , , 



M = 



W 



v/EIm . , 2 /El, . , 

sin k ./ — t • Sin k x 

n n V III n 



(32) 



n— 1 , 3 > 5 • • • 



To further check this, it can be seen by inspection that using equation 
(31) to determine M from the formula. 



M = El 




) 



equation (32) results. 

The shear equation does not converge when > “ , as is the case 

when classical theory is utilized. This fact shall not be derived here. 

Next, look at the equations for y, Q, and M when goes to zero. 

Notice that this can occur in a beam of small for any mode, n ; or for 

any beam at a large n. 

Substituting equation (27b) in (26), (28), and (29), it can be shown 

that 
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y(x,t) = V 



00 
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Certainly these equations look forbidding^ but their derivation involves 
only algebra^ even though long and complicated^ and shall not be indicated 
here, 



Taking the limits of (33)^ (3^)^ and (35) as L'/r’ goes to zero^ one 



finds 



lim 



0 
r ’ 



y(x,t) = V A - % yg ^ k Sin co t sin k^x ,( 36 ) 
I TT ^ ^ 



n— 1 ^3)3 * * • 



lim 



1,-0 



Q(x,t) = ^y - ^ — \ k sin cos k^x , 

r- n 



(37) 



n=l,3,5. .. 



lim 

0 

r'^ 



M(x,t) = 



4tL >/mGA 



IT 



2 — h 

n=l ^3^5* • • 



(^ ■ E ) 



sin t 
G \ 2n 



(^- 1 ) 



sin 



sin k X 
n 



( 38 ) 



The frequencies in the above equations are found by 



lim 

L 



2GAk 



0) 



n 



n 



m 



0 



((^ ^ I) ■ I)) 



from vhich one sees that 
2 



EAk 






n 



GAk 



2n 



m 






n 



In 



m 



(39 a,b) 



One might ask why such complicated equations as (33)^ (3^)^ and (35) 
were used to obtain the limits instead of merely substituting (39 a^b) 
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into (26), (28), and (29). Had this been done, the sin portion 
of the moment equation wo\ild have gone to zero, yielding a faulty 
equation. 

Equations (39 a?h) indicate that is associated with the bending 
of the beam since it contains Young’s modiilus. Similarly, is asso- 
ciated with the shearing action. This holds true for all modes in beams 
with small L/r* and for any beam in the higher modes. One could con- 
sider that the discrepancy in the lower modes is caused by the inter- 
action of shear and bending. There is not sufficient time for the 
modification to occur in the higher modes. 

With this Interpretation, it appears that the response of a very 
short beam is almost entirely due to shear, which is completely acceptable 
intuitively. Only the bending moment shows any effects, then, of rotary 
inertia, since only it has the sin c Op^ t term. 

A zero-length beam probably has no significance physically, but, 
when it is recalled that as n gets large, the nth term goes to the form 
of (36), (3T)^ and (38) for all beams, they become more interesting. 

Notice equation (3T) in particular. This equation can be transformed by 
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n U TT /_ n 
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\ sin k (x'^J— t) 
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( 40 ) 
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This is just the equation for square waves traveling In from the ends. 



The magnitude of each component Is v ^ . The component traveling 

In from x = 0 Is Initially positive, and the one from x = L is initially 
negative. The waves travel at a velocity, 



velocity = 





Further notice that, after each conponent has traveled the length of 

the beam once, a reflected shear wave of magnitude v /mGA results. 

Thus, the transient response of a shaft could be approximated by this 

square wave, and the accuracy Improves as gets small. 

Lastly, it is noticed that equations (36), (37)^ a^nd (38) prove 

the convergence of the basic equations for y, Q, and M, since, as 

indicated above the terms of each of the equations represent the limit 

as n gets large. Thus for the higher modes, y^ and decrease as 

— 2 and decreases as ^ . 

n ^ 



Concept of an Equivalent Beam . An important characteristic of the 
equations for deflection, shear, and bending moment can be realized if 
they are non-dimensionalized. 

Accordingly, define non-dimensional values with primes. Then looking 
at equation (27b), one finds that 




2 



2 2 



n TT 





( 41 ) 



The non-dimensional time is defined such that 



U)’ t’ = 10 t 

n n 
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whereby 



t« = ^ t 

L \/ m 



ik2) 



From the definition of k , 

n 



t ^ 

^ =L 






Continuing, from equations (33)^ (3^)^ and (35); it is seen that 
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From these equations^ it is clear that the non-dimensional values 
of frequency^ deflection, shear, and bending moment are equal for any 
non-dimensional time or position, for any tvo beams if their values of 
z are the same. 

Recalling that 




it is thus necessary that, for the non-dimensional equations of any tvo 

G Li 

beams to be equivalent, the ratios ~ and — , for one beam must equal 
these ratios for the other. 

Look at the ratio of the modified shear modulus to Young *s modulus. 
As was indicated earlier, 

G = k* G* , 



where k* depends on section shape. Evaluation of k* is complex and beyond 
the scope of this work; however, for a first approximation k* could be 
evaluated by the equation. 




o 



(^T) 



where 

b = width of the beam material at the neutral axis, 

C = height of the extreme fiber of a section measured from the 
neutral axis, 

h = height of any point in a section measured from the neutral axis. 
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See reference pj ^or the exact evaluation of Use of equation (^7) 

yields results accurate to about four per cent for a solid circular cross 
section. 

Clearly there are many variations in shape that viould yield a 
k ’ G * 

constant value of — = — . 



Effect of Damping , Even though the equations derived herein are appli- 
cable only in the elastic region^ hysteresis damping does occur. If 
the equatJ.ons for y^ Q, and M are used in the form of equations {26), 

{28), and ( 29 )^ respectively, there is a possibility that many of the 
modes will become additive at some relatively long time. The result woiold 
be inordinately large magnitudes for the three quantities* of interest. 

Applying damping in an analytical way would be difficult to say the 
least. However, its major effect, i.e. that of reducing amplitude, could 
be indicated approximately without loss of generality for the formulas. 

Accordingly let the deflection energy that is lost each cycle be 
two per cent of the energy of that cycle . From this, it can be 
shown that a damping factor for each term would be approximately 

• 4. 4- -0.00l6l'"n'*^ ^,Q^ 

Damping factor = e ^ ) 

See Appendix D for details of this derivation. 

Suggested Design Method . A method of calculating the response of any 
beam is implied by the concept of an equivalent beam. Using equations 
(A2), (45), and (46), modified by damping factors, one could 



-23- 



obtain values of xnaximum shear, deflection and bending moments and the 

Q 

time at which these occur. For a constant — , these maxima could be 

E ’ 

L G 

plotted versus — , . Several plots for various — would result in a 
set of curves that would be completely general. After a few simple 
calculations to convert from the non-dimensional magnitudes to absolute 
magnitudes, one woiild have a good indication of the stresses experienced 
by the shaft of interest. 
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III. RESULTS 



The equations for deflection, shear force and bending moment, 
modified by the appropriate damping factors, were programmed on the 
IBM 7090 computer for a solid steel, circular beam. The results of 
that computation are indicated in Figures II, III, and IV. These plots 
are for unity excitation velocity; thus one need only multiply the 
actual velocity to obtain the magnitudes desired. 

Figure II shows the time history of y, Q and M for a beam of six 
inch radius and thii*ty-six inch length. Notice that, for a -unit velocity 
input, the deflection curve magnitude coiold be considered in units of 
time. The origin of each shear and moment cuive is the horizontal 
line representing the position of the end points on the corresponding 
deflection curve. Since the excitation is symmetrical only one-half of 
the beam is represented. 

Figure III is included to obtain a qualitative picture of the com- 
parison of the responses of two beams of different . The solid line 
is for the beam of Figure II while the dashed line represents the response 
of a beam of six inch radius and sixty inch length. The units of this 
figure have been normalized by dividing time by length and y, Q, and M 
by the product of velocity and length. Since, for both beams, the cross 
section is the same, equations (42), (44), (45), and (46) indicate the 
validity of this normalization. Notice that the magnitudes have not been 
non -dimens ionali zed, although that too could have been done. Nevertheless 
the normalizing method employed provides a satisfactory base for coir 5 )arison. 
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Figure IV shows the variation of the magnitude of the maximum 
deflection, shear, and bending with respect to . Here all units 
have been non-dimensionalized in accordance with equations (kk), (^5)^ 
and ( 46 ). The position of the maximum deflection and bending moment 
was, as would be expected, the center of the beam. The position of the 
maximimi shear varied from x' = 0 to about x' = .125 . 
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Figure IV 



Kon-Dimenslonal Maximum Deflection, Shear 
Force, and Bending Moment 




- 29 - 



IV* DISCUSSION OF RESULTS 



Upon examination of Figure II ^ one's attention is immediately drawn 
to the negative motion that precedes the positive deflection as the 
disturbance travels into the beam. The curves for shear and bending 
moment tend to support this, but no physical reasoning is evident. To 
show how this occurs, look at the solutions for y and y, . By putting 
equations (24) and (25) into (18) and (19)^ it is seen that, relative to 
the end points. 
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(50) 



Equations (^9) s-nd (50) show that the coefficients of each term in 
the series is the same as the coefficient of the corresponding term for y, 
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the net deflection, multiplied by a factor. For simplicity the nth term 
of equations (49) and ( 50 ) could be written, 
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As n goes to infinity, 
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(51) 



(52) 



where (l-) indicates that this factor approaches one. Wow, by examining 

equation ( 26 ), the equation for y, it is seen that the coefficient, of the 

E 

sin term is negative. Furthermore, — is always greater than one, 

e,g. for a solid circular shaft it is approximately 3»45* With these 
facts it can be seen that, for a circular shaft. 
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In other words, this shows that the deflection due to shear is greater 
than the net deflection and that the deflection due to bending is initially 
in the opposite direction] When it is recalled that the higher series of 
frequencies is associated with bending the the lower with shear, it fol- 
lows that the bending effects should travel into the beam faster than 
the shear effects. The ratio of these speeds will be on the order of 



(ju, 



U). 



in 



E 

G ' 



Utilizing these concepts. Figure V shows, in qualitative terms, how 
the resioltant deflection occurs. 

Fi glare V 

Qualitative Description of Deflection Components 




It can be deduced from this discussion, then^that the negative shear 
that precedes the expected positive values is due to the negative bending, 
and not the impact shear. This accounts for the fact that the negative 
portion keeps pace with the bending action instead of moving at the 
slower speed of the shear wave. Conversely, the rising, positive bending 
moment is due to the expected positive shear. 

By equations ( 53 ) and ( 5 ^)> it can be seen that the energy contained 
in a beam is much greater than is indicated by the deflection alone. Even 
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though the deflections due to shear and bending are of opposite sign^ their 
respective energies are additive. Thus the actual energy in a beam is 
much greater than classical theory predicts. Lack of time precludes 
deeper investigation of this point to gain quantitative information. 

It was noticed that the curves for shear are relatively uneven and 
ragged. Even though the computations went to the forty-nipth mode^ this 
is probably not enough since the series converges only as i ^ and as that 
only in the higher modes. 

Upon examining Figure the lack of similarity in the two responses 

is noteworthy. This confirms that beams of two different length to radius 
of gyration ratios cannot be compared^ even in a normalized form. Fur- 
ther^ it is noted that^ as predicted^ the shear curves of the beam with 
the longer look less like the ideal square wave than those of the 

shorter beam. 

The accuracy of the curves in Figure IV is doubtful,. Because of com- 
puter time limitations the search for the maxima may have been done with 
the time and position increments too large. For this reason also^ the 
range of time examined may not have included the absolute maxima. However^ 
the curves are included to show the general variation of the indicated 
quantities with . 

The design method suggested in the Procedures section is dependent 
upon the availability of sufficient data from parent models. To obtain 
sufficient data to be generally useful, machine corr5)utation will probably 
be necessary. However, programming the equations derived in this report 
would not be too difficult. An example of how to utilize such derived 
curves is indicated in Appendix E. 
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A further refinement of the results of this work could be accomplished 
by incorporating them into elasticity theory to obtain the actual maxi- 
mum stresses experienced in a beam. A calculation of this sort could be 

included in the computer program previously mentioned and a single family 

L G 

of curves of maximum stress versus — , for a range of — would result. 

r ill 

Lastly^ it is recognized that the results described herein should be 
confirmed by experiment before any actual engineering use of this is 
attempted. Such experiments are therefore recommended. 
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V. CONCLUSIONS 



It is felt that the equations derived speak for themselves and require 
no further qualification. The equations show that the deflection due to 
shear predominates in the total deflection and that the deflection due to 
bending is in the opposite direction. The energy contained in a beam is 
greater than is indicated by the net deflection because of this fact. 

The equations imply that the disturbing effects of bending travel through 
the beam faster than those of shear. 

It has been shown that beam response may be placed in a non-dimensional 
form whereby all beams with the same ratio of length to radius of gyration 
and the same ratio of shear modulus to modiilus of elasticity, modified by 
a factor accovinting for section shape, will have identical non-dimensional 
responses . 

This fact forms the basis of a siii5)lified method to obtain design 
calcinations. The method is, however, dependent upon the compilation of 
response data from a number of parent models. 
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VI. RECOMMENDATIONS 



It is recoramended that the equations derived in this work be 
confirmed by experiment. 

The indication that significantly greater energies are contained 
in a beam than is implied by the net deflection should be investigated 
further . 

It is recommended that graphs for the non-dimensional form of the 
response of any beam be obtained. In doing so, care shovild be taken to 
assure that enough terms in the series be included to provide sufficient 
accuracy and to assvtre that enough of a span of time is examined, so that 
certainty of obtaining the maxima resvilts. 

Lastly, it is recommended that the results obtained in this work 
be applied to the theory of elasticity so that prediction of the magnitude, 
time and position of maximum stress co\ild be calciilated. It is recommended 
that the results of this be placed into a family of curves, where each 
curve is for a constant modified shear modulus to Young's modulus ratio, 
and the coordinates are maximum stress and the ratio of length to radius 
of gyration. 
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Appendix A, Application of Boundary Conditions 



Equation (8) specifies that X(0) = X(l) = 0. However this could 
be specified more completely by, 



Xg (0) = Xg (L) = X^ (0) = X^ (L) = 0 



Now let = C, + C^, • , and so on. From this, one can look 

1 Is lb ^ 2 2s ^ 2b’ ’ 

at the boundary conditions for the deflection due to bending sepa,rately. Thus, 



X^(?^) = sin kx + cos kx + sirh 'kx + cosh^ikx , (a) 



subject to boimdary conditions 



Xj^ (6) = X^ (L) = 0, 






Sx 






x=0 



3x 



= 0 



x=L 



(b) 



Carrying out the Indicated operations in (a) and (b) above, 
\ (°) ■ =2b " =l.b ■ ° ' 



3x^ 



= -=2b * <=l,b ■ ° 



x=0 



Then 



Therefore C„ = Ci , = 0. 
2b 4b 



(L) = sin kL + sirih '.kL = 0 , 
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Bx 



= - C-, sin kL + C_, sin hkL = 0. 
lb 3b 



x=L 



Since sin hkL can never be zero, the above can only be true if C_, = 0 and 



3b 



kL = nn , n = 1, 2, 3> •••• 

Thus, 

X^(x) = sin 21 X. 



(c) 



But the constants for X , the deflection due to shear, have not yet been 

s 

satisfied. Thus, 



Xg (x) = sin kx + cos kx + sinhlikx + cosh’ikx , (d) 



and 






X^ (L) = sin kL + cos kL + sirh Mi + cosh M = 0 . 



"2s 



3s 



"4s 



But from ( c ) k = ^ , 

^ L , thus. 



Xg (L) = C 2 g + C^g sinh nrr + cosh jiTT 



= C 2 g (l - cosh nrr) + C^g siuh nn = 0 . 



But this is only true if = C_ = 0. Thus, X = (C^ + C^, ) sin 

ciS JS J.S -LD h 



= C sin 



nrrx 



Since the excitation of the beam motion is syinmetrical, the slope of the 
beam at x = ^ must be zero. This precludes any of the even numbered modes, and 



X^ = sin n = 1,3,5,.... 



(e) 
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Appendix B 



More Convenient Form for Frequencies 



Starting with equation (2j), multiply numerator and denominator by, 
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For the case under study, a solid circular beam, 
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Now multiply numerator and denominator of the last term under the radical 
El 

"to obtain 
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2 

Thus, by putting this into the equation for , and factoring. 




(b) 



-hi- 



Appendix C 



Proof of Convergance of Bending Moment Equation for Large 
Length to Radius Ratio . 



From equations ( 29 ) and (30 a,h) , it can be seen that, letting 
be the nth term in the infinite series of M, 
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Look at the expression from the above equation, 
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Then 




With this and equation (27b), since never goes to infinity, 

■^n 




Remembering that. 




it can be seen that (a) goes to the form ^ as goes to infinity. 



Clearly L' Hospital's nile applies here, but, in order to keep the 
application as simple as possible, recall that this calculation is only 
to show that this limit goes to zero. As the numerator is differentiated 
in accordance with L 'Hospital’s rule, it will always have zero as a limit 
as goes to infinity. Thus the case will be proved sin5>ly by showing 

that the limit of some derivative of the denominator becomes other than 
zero. 

To proceed, then. 
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Since we are not interested in the actual magnitude of the limit, cancel 
out factors which are constant with respect to . Thus we get 




Therefore the sin cu^ t term goes to zero in the limit as — , goes 

£:n r 

to infinity. 

By similar reasoning, simply by multiplying J by ^ in the sin 

term. 
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Appendix D, 



Derivation of Damping Factor 



A beam experiencing transverse motion could be interpreted as an 
extremely complex spring-mass system. Since the energy of any spring- 
mass system is proportional to the square of its deflection, one could 
write 

El =< 7 ^ , 

where 

= Energy of the ith cycle, 
y^ = Deflection of the ith cycle. 

If the energy lost per cycle is two per cent of the energy of that 
cycle, then 
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(a) 



From elementary vibration theory the damping factor is of the form, 

-00 Bt 

Damping factor = e ^ , (b) 

where 

B = fraction of critical damping, 
and the frequency, modified by damping is 

00, = 00 v/l - B^ . (c) 

d n V 

From (c) the time for one cycle, or the period, is 
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Clearly, then, from (b) and (d) , 

-2ttB 

^1+1 /l-B^ 

= e = e'' 



If B is small compared to 1 , 




(e) 



Conibinlng (a) and (e) one obtains 



'i+1 



= 1 - e 



4ttB 



= .02 



Solving this for B one obtains 

B = 0.00161 . (f) 

Notice that the effect of B on the natural frequency in (c) above 
is small. Therefore this effect shall be neglected. 

If each term in the equations for y, y , and y is multiplied by 
its appropriate damping factor and the equations for Q and M are derived 
again, it can be seen that, since B < < 1 , the same equations result 
approximately except that each term is multiplied by the appropriate 
damping factor. 
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Appendix E. 



Example Problem 



Assiame we have a solid steel circular shaft l8 inches in diameter. 
What would be the maximum length that could be allowed between supports 
if the beam is not to deflect more than .1 inch when subjected to a shock 
velocity of 35 inches per second? 

Let steel have the characteristics. 



p = 7-33 X 10"^ 

£ 

G' = 11.2 X 10 



lb. -sec. 

ir 



in. 



psi 



It can be seen that, by equation (hj), 
2 2 

^ hr' ^ R X R ^ 3 

’ / hdA 2h /R^-h^ dh ^ 

Jo Jo 



Thus 

G = ^ G' = 8.45 X 10^ psi. 

By equation (44), 

IF = -1 /8.45 X 10^ ^ ^ 

VP ^ ^ V 7.33 X lO"^ ^ 

For a circular cross section, = 2 ' 

Tabulate y' and — , for several values of L, as indicated below: 

L (inches) y' (^^ p (j^) 

96 3.1 21.3 

120 2.56 26.7 

2.13 
-47- 
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Figure Vl 

Determination of Maximum Length of Beam 
According to Exan5)le Problem 
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Plotting these values in Figure VI,' , which is just a copy of Figure IV, 
it can he seen that the required L/r' is 30' 3» Then 

L = X r’ = 30*3 X 4.5 = 136 inches. 

The maximiim bending moment for a beam with this L/r' is obtained 
from Figure VI and equation h6, i.e., 



M = M' (vL \/mGA ) = M' (vLA ) 

= 35 X 136 X TT X 9^ /t'33 X lo"^ x 8.45 x 10^ 

= 419,000 in. -lbs. 

Using the formula for maximum bending stress, 
c MC 



we find 



S = = 723 psi . 

Similarly, from equation (45) and Figure VI, 



Q = Q' (v /mGA ) = .73 x 35 x tt x 9^ / t -33 x 10“^ x 8.45 x 10^ = 2010 lb. 



Now, the maximum shear stress is 



S 



s 



i. Q 

k’ A 




2010 

'^2 

TT X 9 



= 10.52 psi . 



Clearly the deflection governs in this case, but if the allowable 
deflection had been greater, the acceptable length would have increased. 
This would cause a rise in bending stress. A point would thvis be reached 
where the stresses govern the length. 
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